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We have extended the procedure to find the exact solution of the Riemann problem 
in relativistic hydrodynamics to a particular case of relativistic magnetohydrodynamics 
in which the magnetic field of the initial states is tangential to the discontinuity and 
orthogonal to the flow velocity. The wave pattern produced after the break up of the 
initial discontinuity is analogous to the non magnetic case and we show that the prob- 
lem can be understood as a purely relativistic hydrodynamical problem with a modified 
equation of state. The new degree of freedom introduced by the non-zero component 
of the magnetic field results in interesting effects consisting in the change of the wave 
patterns for given initial thermodynamical states, in a similar way to the effects arising 
from the introduction of tangential velocities. Secondly, when the magnetic field domi- 
nates the thermodynamical pressure and energy, the wave speeds approach the speed of 
light leading to fast shocks and fast and arbitrarily thin rarefaction waves. Our approach 
is the first non-trivial exact solution of a Riemann problem in relativistic magnetohy- 
drodynamics and it can also be of great interest to test numerical codes against known 
analytical or exact solutionsf . 



1. Introduction 

The decay of a discontinuity separating two constant initial states (Riemann problem) 
has played a very important role in the development of numerical codes for classical 
(Newtonian) hydrodynamics after the pioneering work of Godunov (1959). Nowadays, 
most modern high-resolution shock-capturing methods (LeVeque 1992) are based on the 
exact or approximate solution of Riemann problems between adjacent numerical cells and 
the development of efficient Riemann solvers has become a research field in numerical 
analysis in its own (see, e.g., Toro 1997). The success of high-resolution shock-capturing 
methods in many areas of computational fluid dynamics has triggered their extension to 
classical magnetohydrodynamics (e.g., Brio & Wu 1988; for an up-to-date discussion of 
the issue, see Balsara 2004). 

As in other fields in physics, during the last two decades astrophysics, relativity and 

f The code computing the exact solution is available under request to the authors. Users of 
the code can give credit by mentioning the source and citing this paper. 
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cosmology have become computational sciences. Modeling and understanding fluid dy- 
namics in astrophysical scenarios is now a key part in research projects involving su- 
pernovae, relativistic jets, neutron star instabilities, or accretion onto compact objects, 
who share a common distinctive feature: either special or general relativity effects are 
relevant. With this motivation, Ricmann solvers have been introduced in numerical rel- 
ativistic hydrodynamics since the beginning of the nineties (Marti, Ibaficz & Miralles 
1991). Presently, the use of high-resolution shock-capturing methods based on Ricmann 
solvers is considered the best strategy to solve the equations of relativistic hydrodynamics 
which has caused the rapid development of Ricmann solvers for both special and general 
relativistic hydrodynamics (see, e.g., Marti & Miiller 2003, Font 2003). 

The finding by Marti & Miiller (1994) of the analytical solution for initial states where 
the flow is normal to the initial discontinuity boosted the efforts to develop exact Rie- 
mann solvers for relativistic hydrodynamics. Pons, Marti & Miiller (2000) extended the 
domain of solutions to problems with arbitrary initial velocities. Later, Rezzolla & Zan- 
otti (2001), for purely normal flow, and Rezzolla, Zanotti & Pons (2003), for the general 
case, proposed a new procedure to find the solution of the Riemann problem that uses 
the relativistically invariant relative velocity between the unperturbed initial states. How- 
ever, up to date, no analytical or exact solution of the equations of the Riemann problem 
in relativistic magnetohydrodynamics has been derived. 

The equations of both classical and relativistic magnetohydrodynamics form a non- 
strictly hyperbolic system. A consequence of the non-strict hyperbolicity of the MHD 
are the degeneracies in the wave speeds (that lead, e.g., to compound waves, admisiblc 
solutions of the planar MHD that involve intermediate shocks), which must be han- 
dled analytically with care and hinder the development of both exact and approximate 
Riemann solvers using the characteristic information for the MHD equations. In the 
relativistic case the difficulties in the development of such solvers are increased by the 
higher non-linearity of the RMHD system of equations. The characteristic structure of the 
equations of RMHD is analyzed in the books by Lichnerowicz (1967) and Anile (1989). 
Approximate Riemann solvers using the characteristic information have been developed 
by Romero ct al. (1996), for the same particular magnetohydrodynamic configuration 
considered here, and by Balsara (2001), Komissarov (1999) and Koldoba, Kuznetsov & 
Ustyugova (2002), for the general case. A number of analytical solutions involving only 
shocks, rarefactions and Alfvcn waves (Komissarov 1999, Komissarov 2003) have also 
been recently derived. 

In this paper we describe the solution of the Riemann problem for the particular case 
in which the flow speed has two non-vanishing components and the magnetic field is 
orthogonal to them. Besides this, we force the flow to have dependence on one spatial 
coordinate taken along one of the two non-vanishing velocity components. For this par- 
ticular setup, the Ricmann structure degenerates to only three waves making the solution 
attainable. The solution reveals interesting and distinct properties of RMHD and could 
serve as a guide in the way to the general RMHD Riemann solution. 

In this paper we closely follow the structure and notation used in Pons, Marti & 
Miiller (2000; hereafter PMM). The paper is organized as follows. Section 2 collects the 
relevant equations. Sections 3 and 4 describe, respectively, the flow across rarefactions 
and shocks setting the ingredients for the Riemann solution, which is discussed in Sect. 5. 
Conclusions are gathered in Sect. 6. 
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2. Equations 

Let JP, T^ v and F*^ u (fi,v = 0,1,2,3) be the components of the density current, 
the energy-momentum tensor and the Maxwell dual tensor of an ideal magneto- fluid, 
respectively 

J»=pu» (2.1) 
= phuV + rp v p - WW (2.2) 
F*^ = u^W - u v W (2.3) 

where p is the proper rest-mass density, h = 1 + e + p/ p + b 2 / p is the specific enthalpy 
including the contribution from the magnetic field (b 2 stands for Wb^), e is the specific 
internal energy, p the thermal pressure, p = p + b 2 /2 the total pressure, and rf 111 = 
diag(—l, 1,1,1) the Minkowski metric in Cartesian coordinates. Throughout the paper 
we use units in which the speed of light is c = 1. 

The four-vectors representing the fluid velocity and the magnetic field in the fluid rest 
frame, u M and W, satisfy the conditions u^u^ = — 1 and u^b^ = 0, and there is an 
equation of state p = p(p, e) that closes the system. All the discussion will be valid for a 
general equation of state but results will be shown for an ideal gas, for which p = (7— l)pe, 
where 7 is the adiabatic exponent. 

The equations of ideal RMHD correspond to the conservation of rest mass and energy- 
momentum, and the Maxwell equations. In flat space-time and Cartesian coordinates, 
these equations read: 

= (2.4) 
= (2.5) 

F*^ = (2.6) 

We consider a particular case in which the flow speed has two components and the 
magnetic field is orthogonal to them. Besides this, we force the flow to have dependence 
on one spatial coordinate (x) taken along one of the two non-vanishing velocity com- 
ponents. Specifically we set u M = W(l, v x , 0, v z ), W = (0,0,6,0), where W is the flow 
Lorentz factor. With these restrictions, the above system can be written as a system of 
conservation laws 

where 

U = (D, S x , S z ,f, B) T (2.8) 
is the state vector of conserved quantities and 

F = {Dv x , S x v x + p, S z v x , S x ,Bv x ) T (2.9) 
is the corresponding vector of fluxes, with 

D = pW (2.10) 

5* = phW 2 v\ (i = x,z) (2.11) 

and 

T = pfiW 2 -p (2.12) 
being the rest-mass, momentum and total energy densities, and 



B = bW 



(2.13) 



4 



R. Romero et al. 



the y-component of the magnetic field as measured in the laboratory frame. Hence, 
according to these equations, the particular initial configuration chosen together with 
the imposed symmetry prevent the generation of new components of the velocity and 
magnetic field. 

It is worth noting that for the particular configuration chosen, the term (6 M 6 I/ ), M ap- 
pearing in the equation of conservation of the stress-energy tensor, vanishes and the 
RMHD equations reduce to the purely hydrodynamical case with the only contributions 
from the magnetic field appearing in the pressure and specific enthalpy, and an addi- 
tional continuity equation for the evolution of the transversal magnetic field. This fact is 
considered in the Appendix where we explore the possibility of including the magnetic 
effects of the present configuration in the definition of the equation of state. 

According to the previous discussion, the magnetized flow under consideration falls in 
one of the two degeneracies of the RMHD system (Degeneracy I; Komissarov 1999), for 
which a description in terms of just three characteristic waves (namely the entropy wave 
and the two fast magnetosonic waves) is adequate. Turning now towards the solution of 
the Riemann problem in this particular case, the discontinuity in the initial states breaks 
down into a couple of left and right propagating rarefaction waves (self-similar contin- 
uous flows) and/or shocks and a central tangential discontinuity across which the total 
pressure, p, is constant. Both thermal and total pressure increase at fast magnetosonic 
shocks. Hence, we would use the comparison of the total pressure at the two limiting 
states to select between shocks and rarefaction waves. 



3. Flow across rarefactions 

Rarefaction waves are self-similar solutions of the equations that depend on x and t 
only through the combination £ = x/t. By imposing such a dependence in system (2.7) 
above we obtain the following set of equations 

(v x - 0^ + P(l + v*W*(v* 0)^ + pW*v'(v* - 0^ = (3.1) 

W*ph(v* - 0^ + 6(1 - ^)| + (1 - «-0| = (3.2) 
W 2 ph(v* - 0^ - Vbt^ - ^| = (3.3) 

rJu x fln z clh 

6(1 + v x W 2 {v x - £))— + bW 2 v z (v x -£)_ + („*- £)— = 0, (3.5) 

similar to the system obtained in PMM. The quantity c s is the sound speed, defined by 

(3.6) 



Cs = 



'1 dp 
h dp 



where s is the specific entropy and h = 1 + e + p/p, the specific enthalpy. 

Non-trivial similarity solutions exist only if the determinant of system (3.1)-(3.5) van- 
ish. This leads to the condition 



± _ - c 2 ) ± aV(l - g) [1 - gg - - gj] 

4 - 1 - [6J) 
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where to 2 = c 2 +v\ — (? s v\ and v\ = b 2 /ph is the Alfven velocity. The plus and minus signs 
correspond to rarefaction waves propagating to the left 1Z^- and right 1Z~> , respectively. 
Note that the values of £ reduce to those obtained in PMM by replacing lo by c s (i.e., 
6 = 0). 

From system (3.1)-(3.5), after some algebraic manipulations, we get: 

W 2 ph(v* ~ 0^ + 6(1 " + (1 - ^)| = (3.8) 

^-hc 2 ^ (3 9) 

kWv z = constant (3.10) 
b 

- = constant. (3.11) 
P 

Now, using (3.9) and (3.11) to eliminate the differentials of b and p, and defining 
B = b/p, the ODE (3.8) can be rewritten as 

<w = ( 1 + €f) (i-fo*) , 3 12) 

dp phW 2 (£-«*)' 

in complete analogy with Eq. (3.20) in Rezzolla, Zanotti & Pons (2003) in the case 6 = 0. 
The analogy can also be extended to the solution procedure. If we define A = hWv z 
then, from Eq. (3.10), 

( v ^=A 2 ( 1 r {v T ), (3.13) 



h 2 +A 2 

which allows to eliminate the dependence on v z in (3.7). Now, from the definition of the 
Lorentz factor, it can be derived 

h 2 4- A 2 

W 2 = . h +A (3.14) 

h 2 (l -(v*) 2 ) 

and, after some algebra, 



Jh 2 +A 2 (l - lo 2 ) 

\ ^ = ±-¥ =r- (3.15) 

Z~v x hw K ' 

(where lo is the positive root of lo 2 ). Finally, substituting these last two expressions in 
(3.12), we get 



dv* , (l + ^h 2 + A 2 (l-Lo 2 ) dp 

_± — ^ — ■ ( 3 - 16 ) 



l-(v x ) 2 h 2 + A 2 (*>' 

The left hand side of this expression can be integrated analytically and the right hand side 
involves only thermodynamical variables and constants. Considering that in a Riemann 
problem the state ahead of the rarefaction wave is known, the integration of (3.16) allows 
one to connect the states ahead (a) and behind (6) the rarefaction wave. The normal 
velocity behind the rarefaction wave can be directly obtained as 



vl = tanhC, 



(3.17) 
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where 

e= i lo ji±<) ± r<i±MEI*^* ,3, 8) 

2 & V1-^/ J Pa h?+A* ( ' 

The differential of p in the last integral is taken along the adiabats of the equation of 
state. The isentropic character of rarefaction waves fixes the entropy to that of state a, 
s a . Having this in mind, the ODE can be integrated, the solution being only a function 
of pi,- This can be stated in compact form as 

v x b =nUVb). (3.19) 

It is interesting to have an expression for the normal velocity inside the rarefaction 
wave in terms of the total pressure. This expression can be built taking into account 
that, from the definition of p, 

dp = dp + B 2 p dp (3.20) 

and that along an adiabat, 
which combined result in 

Substitution in Eq. (3.18) gives 



(3.21) 



d P = 1 1 + t^t I d P- ( 3 - 22 ) 




Note that the previous expression is identical to the one derived by Rezzolla, Zanotti & 
Pons (2003) in the case of (non-magnetized) relativistic hydrodynamics after removing 
the hats and substituting u by the sound speed. The corresponding compact notation for 
the function giving vf in terms of pb will be 

vt=KUPb)- (3-24) 

Function TZ^±(p) is shown in Fig. 1, for different values of the invariant B, the var- 
ious branches of the curves corresponding to rarefaction waves propagating towards or 
away from a. Rarefaction waves move towards (away from) a, if the pressure inside the 
rarefaction is smaller (larger) than p a . The last assertion also applies for the total pres- 
sure, p. In a Riemann problem the state a is ahead of the wave and only those branches 
corresponding to waves propagating towards a in Fig. 1 must be considered. Moreover, 
one can discriminate between waves propagating towards the left and right by taking 
into account that the initial left (right) state can only be reached by a wave propagating 
towards the left (right). The addition of a transverse magnetic field in the limiting state 
forces the value of the normal velocity within the rarefaction wave to have larger absolute 
values. This effect is a consequence of the fact that the absolute value of the slope of the 
function TZ^±(p) at p a , \TZ^±(p a )\, is an increasing function of B (or b a ). Moreover it can 
be easily proved that 

|ft^'(p ;B-oo)| = —\H a ^(p a ;B = 0)\, (3.25) 

where c sa is the sound speed at state a. Hence the sound speed at state a limits in practice 
the range of values of the normal velocity within the rarefaction wave (see Fig. 1 where 
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-1.0 -0.5 0.0 0.5 1.0 

Figure 1: Loci of states which can be connected with a given state a by means of 
relativistic rarefaction waves propagating to the left (71^) and to the right (7\L^) and 
moving towards or away from a. Solutions for B = 0, 1, 3, 10, correspond to solid, dashed, 
dashed-dotted and dotted lines, respectively. The state a is characterized by p a = 0.6, 
p a = 1.0, and — —0.3. An ideal gas EOS with 7 = 5/3 was assumed. 

the curves corresponding to B = 3 and 10 almost coincide). Another consequence of 
the previous result is that in the extreme case for which c sa tends to the light speed, 
the presence of a transverse magnetic field in the limiting state would have no practical 
effect on the rarefaction wave. The effect of the magnetic field must be combined with 
the one coming from the presence of tangential velocities in state a which operates in the 
opposite direction (see PMM). 



4. Jumps across shocks. 

If T, is a hyper-surface in Minkowski space time across which p, and F*^ v arc 
discontinuous, the Rankinc-Hugoniot conditions are given by (Lichnerowicz 1967, Anile 
1989) 

[p^K = , (4.1) 

[T">„ = , (4.2) 

[F*H n v = (4.3) 
where is the unit normal to S, and where we have used the notation 

[G] = G a - G b , (4.4) 



8 



R. Romero et al. 



G a and G b being the boundary values of G on the two sides of S. 

Considering £ as the hyper-surface in four-dimensional space describing the evolution 
of a shock wave normal to the x-axis, the unitarity of n u allows one to write it as 



n v = W s (V s , 1,0,0), 



(4.5) 



where V s is interpreted as the coordinate velocity of the surface that defines the position 
of the shock wave and W s is the Lorentz factor of the shock, 



W s 



1 



(4.6) 



Equations (4.1) and (4.3) allow one to introduce two invariants across the shock 

j = W s D a (V s - <) - W s D b (V s - <), (4.7) 

/ = W s B a (V s - vt) = W s B b (V s - <), (4.8) 

where B — bW . Quantity j represents the mass flux across the shock and according to 
our definition, j is positive for shocks propagating to the right (the same convention as 
the one used in Marti & Miiller 1994, and PMM). Dividing Eq. (4.7) by Eq. (4.8), we get 
that the quantity B/D (or, equivalently, b/p) is constant across shocks as it was through 
rarefaction waves. 

Next, the Rankinc-Hugoniot conditions (4.1), (4.2) can be written in terms of the 
conserved quantities D, and f , and j as follows 

" 1 " 



v- = - 



W s 



D 



\v\ = — 



~D 



D 

= 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



Now from Eq. (4.11) we have that the quantity hWv z is constant across the shock, as 
it is through rarefactions. 

We note that in deriving equations (4.9)-(4.12) we have made use of the fact that the 
mass flux is nonzero across a shock. The conditions across a tangential discontinuity imply 
continuous total pressure and normal velocity (by setting j = in equations (4.9), (4.10) 
and (4.12)), and an arbitrary jump in the tangential velocity and transverse magnetic 
field. 

Our aim now is to write v£ , the normal flow speed in the post-shock state, as a function 
of the post-shock pressure p b . As a first step, we write as a function of pi, j and V s 
(and the preshock state, a). Given the complete analogy between the jump conditions 
Eqs. (4.7), (4.9)-(4.12), and the corresponding expressions in PMM, we write 



vt = h a W a v x a + 



W s (p b -p a ) 



h a W a + {p b -p a ) 



W s v x a 

j 



+ 



1 



PaW a 



(4.13) 



The dependence on the magnetic field in the pre- and post-shock states is hidden in the 
definitions of p and h. 
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The shock speed V s can be eliminated using the definition of mass flux to obtain 
± _ pjWlvl ± \jWf + pjW*{l = vf) 

where (V~) corresponds to shocks propagating to the right (left). 

Proceeding in the same way as in PMM (i.e., [T^] n l/ {{hu^) a + (hu^b} = 0) to derive 
the Taub's adiabat, we can now obtain 

(4.15) 

i.e, the Lichnerowicz's adiabat (Anile 1989) particularized to our special setup. Figure 2 
represents the function 

C a {kh) = h 2 -hl - [ —A— + — ] (p b -Pa) (4.16) 





for an ideal gas equation of state although the general shape of the curve (positive 
asymptotic branches; negative value for h = 1) is independent of the equation of state. 
The (unique) root at the right of h = 1 defines the thermodynamical post-shock state. 

Equation (4.15) together with the definitions of p and h, the equation of state and 
the constancy of bj p through the shock, allows to write pb as a function of pb and the 
preshock state a. Next, multiplying (4.2) by n M and using the definition of relativistic 
mass flux one obtains 

,2 _ - [P 



r 



hjp 



(4.17) 



Using the positive (negative) root of j 2 for shock waves propagating towards the right 
(left), equation (4.17) allows one to obtain the desired relation between the post-shock 
normal velocity and the post-shock pressure pb- In a compact way the relation reads 

vZ=S^(p b ). (4.18) 

Alternatively, the relation can be written as a function of pb 

v* b =SUh)- (4.19) 

Let us note that the expressions used to build up the function S^±, namely Eqs. (4.13), 
(4.14), (4.15) and (4.17), are formally identical to those corresponding to the pure (i.e., 
non-magnetized) relativistic hydrodynamical case. The difference appears in the defini- 
tion of the function p = p(h,p), which leads to different roots of the function C a (h;pb), 
Eq. (4.16). Function S?±(p) is shown in Fig. 3, for different values of the invariant B, 
the various branches of the curves corresponding to shock waves propagating towards 
or away from a. In order to select the relevant branch of the function S^±(p) the same 
argumentation as in the case of rarefaction waves can be used (see §3). As in the case of 
rarefaction waves, the addition of a transverse magnetic field in the limiting state forces 
the value of the normal velocity in the pre/post shock state to have larger absolute values. 
Again, this effect must be combined with the one coming from the presence of tangential 
velocities in state a (see PMM). 
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-10 1 2 3 4 5 6 

A 

h 

Figure 2: Graphical representation of the function C a (h;fib) (definition: see text) whose 
zeroes for varying post-shock pressures, fib, define the Lichnerowicz's adiabat. The state a 
is characterized by p a — 0.25, p a = 1.0 and — 0.5. Solutions for transversal magnetic 
fields in state a, B a = 0,0.5,1.0,2.0, correspond to solid, dashed, dashed-dotted and 
dotted lines, respectively, fib was chosen to be 1.0. An ideal gas EOS with 7 = 5/3 was 
assumed. 

Once vl is know, can be obtained through 

,,20, 

where we have defined A — h a W a v*- Analogously, bb — Bpb with B = b a /p a . 

5. The solution of the Riemann problem. 

As discussed in § 2, for the particular case under consideration (magnetic field orthog- 
onal to both the fluid velocity and the wave propagation direction), the time evolution 
of a Riemann problem with initial states L (left) and R (right) can be represented as: 

LR -> LW^ L*C R*W^ R (5.1) 

where W and C denote a (fast magnetosonic-) shock or rarefaction, and a contact dis- 
continuity, respectively. The arrows (<— / — ») indicate the direction (left / right) from 
which fluid elements enter the corresponding wave. 

The solution of the Riemann problem consists in finding the intermediate states, L* 




Figure 3: Loci of states which can be connected with a given state a by means of 
relativistic shock waves propagating to the left (<S«_) and to the right (<S_ >) and moving 
towards or away from a. Solutions for B = 0, 1, 3, 10, correspond to solid, dashed, dashed- 
dotted and dotted lines, respectively. The state a is characterized by p a — 0.25, p a — 1.0, 
and — 0.5. An ideal gas EOS with 7 = 5/3 was assumed. 



and as well as the positions of the waves separating the four states (which only 
depend on i, L„, _R* and R). The functions and W«_ allow one to determine the 
functions v^(p) and respectively. The pressure^* and the flow velocity in the 

intermediate states are then given by the condition 

«L(P.) = «L(P.) = <■ (5-2) 
The functions Vg^{p) are defined by 

K s (p) tfp<p s 
S s (p) tfp>p s 



sM = \ ~ s :< ' (5-3) 



where TZ s (p) (S s (p)) denotes the family of all states which can be connected through a 
rarefaction (shock) with a given state S (L, R) ahead of the wave. Once p* and v% have 
been obtained the remaining quantities can be computed. 

Figure 4 shows the solution of a particular Ricmann problem for different values of 
the magnetic field 6 = 0, 1.0, 2.0, 4.0 in the initial states. The crossing point of any two 
lines gives the pressure and the normal velocity in the intermediate states. 

It must be noted that the resolution of the Riemann problem under consideration can 
be formally done in the same way as the pure (relativistic) hydrodynamical problem with 
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Figure 4: Graphical solution in the (j>, i^-planc of the rclativistic Riemann problem 
with initial data p L = 1.0, pl = 1-0, v x L = 0.0; pr = 0.1, pjj = 0.125 and = 0.0 
for different values of the transversal magnetic field B = 0,1.0,2.0,4.0, in the left and 
right states represented by solid, dashed, dashed-dotted and dotted lines, respectively. 
Diamonds indicate the initial states. An ideal gas EOS with 7 = 1.4 was assumed. 
The crossing point of any two lines gives the pressure and the normal velocity in the 
intermediate states. Iq gives the solution for vanishing magnetic field. 

some modifications (see the Appendix). First of all, p and h have to be replaced by p 
and h. Secondly, the sound speed in the integrand of Eq. (3.18) has to be replaced by 
lo. Finally, the equation of state that provides the rest-mass density as a function of the 
pressure and the enthalpy, p(h,p), has to be modified to include the contributions from 
the magnetic field. Given the parallelism between the present particular (relativistic) 
magnetohydrodynamical case and the purely hydrodynamical one, the effects concerning 
the smooth transition from one wave pattern to another when the tangential velocities in 
the initial states are changed (Rczzolla & Zanotti 2002, Rczzolla, Zanotti & Pons 2003) 
will extend to the present case for fixed initial values of the magnetic field. Hence we 
concentrate in the effects on the solution induced by varying the initial magnetic fields. 
Figure 5 shows the solution of a Riemann problem with a) vanishing magnetic field, and 
b) 6_r = 0.8. Whereas in the purely hydrodynamical case the Riemann solution gives rise 
to a left propagating rarefaction wave and right propagating contact and shock waves, the 
case with non vanishing magnetic field leads to a couple of shock waves and a left propa- 
gating contact discontinuity. The reason for this qualitative change (rarefaction/shock to 
shock/shock) can be found in the increase of the total pressure in the initial right state 
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of the magnetized case. This increases the total pressure in the intermediate states at 
the two sides of the contact discontinuity. When this pressure becomes larger than that 
at the initial left state then a left propagating shock instead a rarefaction is produced. 

Also noticeable from Fig. 5 is the increase of velocity of the shock propagating towards 
the right in the magnetized case. The increase of the velocity of propagation of fast 
waves for increasing magnetic fields (approaching the light speed in the fluid rest frame 
for strong enough fields, much larger than equipartition) is well-known in RMHD (e.g., 
Anile 1989). In the particular magnetic problem under consideration it leads to the 
increase of velocity of propagation of rarefaction heads and tails, and of shocks and can 
be understood as follows. Equation (3.7) gives the propagation speed of the head/tail of 
a right/left propagating rarefaction wave on a given state. Taking v x ,v = in Eq. (3.7), 
one gets = lo and u — > 1 for large enough b. One should remember that for purely 
hydrodynamical (relativistic) flows, | has the sound speed as limiting value. A similar 
result holds for shocks. For preshock states at rest, \V^\ — > 1 as the magnetic field in the 
preshock state is increased. This can be seen by remembering that the shocks that appear 
in our configurations are super-magnetosonic, and that the fast magnetosonic speed in 
the preshock state (uj a ) tends to light speed when b a — > oo. 

Finally, let us note that the drift towards solutions involving only discontinuous waves 
(shock waves and rarefactions of negligible width) for increasing magnetic fields as con- 
cluded in the previous paragraph, is consistent with the fact that in the limit of strong 
magnetization, the equations of RMHD reduce to the equations of force-free degenerate 
electrodynamics, whose Riemann problem only involves (linearly degenerate) discontin- 
uous waves (Komissarov 2002). 

6. Summary and conclusions 

We have obtained an exact solution of the Riemann problem for multidimensional 
relativistic magnetohydrodynamics in the particular case in which the magnetic field is 
normal to the fluid velocity. In this particular problem, the complex 7-wave pattern of 
RMHD is reduced to two fast magnetosonic waves and a contact discontinuity, which 
allows to use the same procedure as in the non magnetic case. Alternatively, we have 
shown that the problem can be understood as a purely RHD situation with a modified 
equation of state (see Appendix A for details). 

Two interesting features arise from our results. First, for fixed initial thermodynam- 
ical states, it is possible to change continuously from one wave pattern to another 
(shock/shock, shock/rarefaction, rarefaction/rarefaction) analogously to what happens 
when tangential velocities are introduced (Rezzolla & Zanotti 2002). Secondly, we re- 
cover the result for RMHD flows with general magnetic field configurations stablishing 
the tendency of the wave speeds to the light speed when the magnetic field dominates the 
thermodynamical pressure and energy. For our particular configuration of the magnetic 
field, this results in fast moving shock waves and rarefaction waves in which the distance 
between the head and the tail is progressively reduced as the magnetic field increases. The 
drift towards solutions involving only discontinuous waves (shock waves and rarefactions 
of negligible width) for increasing magnetic fields is consistent with the fact that in the 
limit of strong magnetization, the equations of RMHD reduce to the equations of force- 
free degenerate electrodynamics, whose Riemann problem only involves discontinuous 
waves (Komissarov 2002). 

In addition to the theoretical interest of our results by their own, having obtained 
an exact solution of the RMHD Riemann problem is relevant for the development of 
numerical codes. Up to now, in order to test the various algorithms and approximate 
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Figure 5: Solution of the Ricmann problem with initial data pl — 1.0, pl = 1-0, 
v x L = 0.0; p R = 0.1, p R = 0.125 and v x R = -0.5, at time t = 0.4. An ideal gas EOS with 
7 = 5/3 was assumed. Left panel: purely hydrodynamical problem. Right panel: with 
bR = 1.0. Note that p is constant through the contact discontinuity whereas the thermal 
pressure, p, is not. 



Riemann solvers developed for numerical applications, one could only increase the spatial 
resolution and hope that the numerical solution converged to the physical one. Having 
an exact solution to compare with, even if it is just a particular case, allows for a more 
rigorous testing and error estimation. Last but not least, from a pedagogical point of 
view, it is more convenient to start understanding and solving a simpler case before 
attempting the solution of the full problem, which is the next natural extension of this 
work. 
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Appendix A. A hydrodynamical approach 

Equations (2.7) are identical to those for the purely hydrodynamical case by replacing 

b 2 

V — > P = P+y (Al) 
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b 2 

h — >h = h+— , (A 2) 

P 

indicating that a description of the present particular RMHD problem based on a purely 
hydrodynamical approach with a different equation of state may be possible. In this Ap- 
pendix, we explore such a possibility, first suggested by Romero et al. (1996). The key 
point is to eliminate the magnetic field from the equations by building up a thcrmody- 
namically consistent equation of state including the effects of the magnetic field. 
It follows from eqs. (2.7) that 

^ = 0. 

where D /Dt stands for the standard convective derivative, implying that the evolution 
of the fluid elements is along states keeping b/p — constant. Then for a particular fluid 
element, p and h can be written as 

p = p + B 2 p 2 /2 (A 4) 

h = l + e + ^+B 2 p/2, (A5) 

where B is a constant. 

Consistency of the two previous expressions with the definition of h is fulfilled by 
defining a new specific internal energy, i, 

e = e + B 2 p/2. (A 6) 

The evolution of the fluid elements in a perfect fluid is adiabatic. Hence now we look 
for the adiabats of the new equation of state, p — p(p, i). The fact that the evolution of 
the fluid elements keeps B = constant draw us to consider that B is constant along the 
adiabats of the new equation of state, s — constant. We shall use this fact to look for the 
desired relation between the entropies of the two equations of state, s and s. To do this, 
we differentiate eq. (A 6) along transformations keeping s = constant. We get 

di = de + B 2 dp/2. (A 7) 

On the other hand, according to the first law of thermodynamics, for an adiabatic trans- 
formation, 

d£=^dp. (A 8) 

Substitution of p in the previous expression gives 

de = ^dp + B 2 dp/2. (A 9) 

Finally, comparison with eq. (A 7) leads to 

de = ^rdp, (A 10) 

P 

which is formally identical to the variation of internal energy in an adiabatic transfor- 
mation s = constant. Taking into account that the differentials were taken along the 
adiabats of the new equation of state, the conclusion is that the entropy in the new 
equation of state must be a function of the entropy in the original equation of state only. 
Now the sound speed of the new equation of state, 
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can be derived. The substitution of p following eq. (A 4) and the equivalence of the 
adiabats of the two equations of state leads to 



where h and c s stand for the enthalpy and the sound speed of the original equation of 
state, respectively. Finally, a bit of algebra allows us to write, 



v\)cl+v% (A 13) 

where va is the Alfven speed for our particular case in which only one component of the 



magnetic field is non-zero, va — b/y ph = B, - — I . Note that c s coincides with 

V V h + B 2 pJ 

the quantity uj defined in § 2. 

It is interesting to note that although the original equation of state could have a sound 
speed significantly smaller than light speed (e.g., < 1/V3, for an ultra- relativistic non- 
degenerate ideal gas), the sound speed of the new equation of state (that represents the 
true propagation speed of perturbations in our magnetized fluid) approaches the light 
speed for large enough values of B (or b). Finally, notice that the convexity of the EOS 



d 2 p 

is ensured, since B 2 is a positive defined quantity and therefore —-, - 

ap z 



> 0. 
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